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We investigate the cjp inflationary model in the Bianchi-I spacetime using the effective 
spacetime description of loop quantum cosmology to understand the issues of the resolution 
of initial singularity, isotropization, effect of anisotropics on the amount of inflation, and 
the phase space attractors in the presence of non-perturbative quantum gravitational modi- 
fications. A comparative analysis with the classical theory by including more general initial 
conditions than the ones previously considered in the latter is also performed. We show 

^^ ' that, in general, the classical singularity is replaced by a bounce of the mean scale factor in 

CN , loop quantum cosmology. Due to the underlying quantum geometric effects, the energy den- 

^.^' sity of the inflaton and the anisotropic shear remain bounded throughout the non-singular 

^ C^ , evolution. Starting from arbitrary anisotropic initial conditions, a loop quantum universe 

isotropizes either before or soon after the onset of slow-roll inflation. We find a double at- 
tractor behavior in the phase space dynamics of loop quantum cosmology, similar to the 

^^ ' one in classical theory, but with some additional subtle features. Quantum modifications to 

the dynamical equations are such that, unlike the classical theory, the amount of inflation 
does not monotonically depend on the initial anisotropy in loop quantum cosmology. Our 
results suggest that a viable non-singular inflationary model can be constructed from highly 
anisotropic initial conditions in the Planck regime. 
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> , I. INTRODUCTION 
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\^ • The inflationary paradigm provides an excellent explanation of the homogeneity and the flatness 

C^^ ■ problems of the standard model of cosmology and the origin of the observable large scale structure 

Tij- I in our universe. However, various questions about the physics of the pre-inflationary stage of 

^^ ' the universe remain unanswered. In particular, in general relativity (GR), inflationary spacetimes 

have been shown to be past incomplete [l|, and little is known about the initial conditions and 

the geometry of the pre-inflationary epoch. Both of these issues need to be carefully addressed 

in any complete model of inflation. Unlike the issue of the past singularity, which can only be 

r> ' faithfully addressed using a quantum theory of gravity, the issue of the role of pre-inflationary 

C^ . geometry on the onset and predictions of inflation has been previously studied using GR (and 

also in modified theories of GR). In this setting, various investigations have been carried out to 

understand the onset of inflation in the presence of anisotropic shear in the pre-inflationary epoch, 

often by assuming a homogeneous and anisotropic patch of spacetime using Bianchi models ^- 



cn 



251]. Many results point towards isotropization and genericness of inflation starting from initial 
anisotropic conditions, however counter-examples have also been found (see for eg. 2ll. |24|. |25||). 
Since the anisotropic shear diverges in classical gravity close to the initial singularity, it has been 
argued that a consistent treatment of anisotropies in the pre-inflationary epoch is incomplete 
without inputs from quantum gravity [9|, llJ|. Due to these reasons, it is important to study the 



pre-inflationary dynamics in a quantum gravitational anisotropic spacetime to understand the way 
quantum gravity effects resolve the classical singularity, and affect the isotropization and the onset 
of inflation in presence of anisotropies. 
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The goal of this article is to address these issues in loop quantum cosmology (LQC) [26|J, a 
framework to quantize homogeneous spacetimes based on loop quantum gravity (LQG). A key 
prediction of LQC is that the backward evolution of a universe, starting from classical initial data, 
by the quantum Hamiltonian constraint leads to a non-singular bounce 27h29||. The resolution of 
singularity in LQC is a result of the underlying discreteness of the quantum geometry at the Planck 
scale, and its occurrence has been confirmed in a variety of isotropic and anisotropic models (see 
for eg. 3CH43|| ) . including the (f)^ inflationary model in the isotropic setting 4J]. A powerful tool 
in understanding the detailed physics of these models is the effective spacetime description of the 
underlying quantum dynamics 45|, |46|], which using sophisticated numerical simulations has been 
tested to the Planck curvature scale for various matter models for both isotropic 27H3d. |34| . ISg . |3a) 
as well as anisotropic spacetimes 471, |48|] (see Refs.[49!, l50| for a review of numerical methods in 
LQC). The effective spacetime description is derived using an effective Hamiltonian obtained via 
geometric methods of quantum mechanics [SlJ , by assuming states which become semi-classical at 
late times, describing a macroscopic universe in the classical regime. It is found that GR is an 
excellent approximation to quantum dynamics at curvature scales small compared to the Planck 
value, however departures between the two theories become significant in the Planck regime which is 
described by the effective dynamics to an excellent accuracy. Using effective Hamiltonian method, 



bounds on expansion and shear scalars have been derived [52l.l53l|. and generic results on the absence 
of strong singularities in isotropic 5J, |55(] and Bianchi-I spacetime [56] have been obtained. 

In this manuscript, we analyze the dynamics of (jP' inflation model in the Bianchi-I spacetime 
in the effective spacetime description of LQC. We first show that unlike in the classical Bianchi-I 
model, the past evolution in LQC is non-singular. Singularity resolution is confirmed for a large set 
of initial conditions, in confirmation with the generic result on resolution of singularities in Bianchi- 
I spacetime in LQC [5a|. In contrast to the isotropic inflationary models in LQC where there is a 
single bounce of the isotropic scale factor 57h64|| always occurring at the maximum allowed value 
of the energy density pmax = 0.41 /Opi, the resolution of singularity in Bianchi-I spacetime occurs 
via non-singular bounces of the three directional scale factors (01,02,03) occurring at a range of 
values of energy density and anisotropic shear scalar a"^, bounded by pmax (same as in the isotropic 
model) and cj^ax = 11.57//|j. In the isotropic LQC, p = Pmax at the bounce constrains the value 
of the inflaton velocity at the bounce for a given initial value of the inflaton field. On the other 
hand, in Bianchi-I spacetime in LQC, there is no such restriction. The bounce of the mean scale 
factor in Bianchi-I spacetime (defined as o = (010203)^'^ can take place even if the energy density 
is not at its maximum value at the bounce. In this way, presence of non-vanishing shear leads to 
sets of initial conditions which were not allowed in the isotropic spacetime. 

The availability of a non-singular inflationary model in the presence of anisotropics sets the stage 
to answer some important questions which so far could only be addressed with in the limitations 
of the classical theory. These include, the way slow-roll inflation in 0^ model begins, starting from 
highly anisotropic geometries; whether slow-roll inflation is an attractor of trajectories starting 
from arbitrary anisotropic conditions in the Planck regime; and the way non-zero anisotropic shear 
affects the number of e-foldings. Though these issues have been addressed to some extent earlier 
in the classical theory, for completeness and comparative analysis, we analyze both the classical 
and the effective dynamics of LQC for the (j)^ model in Bianchi-I spacetime. In the process, apart 
from answering the above questions in the non-singular anisotropic inflationary model in LQC, we 
gain new insights on some of these issues in the classical theory. Below, we summarize the main 
results on the study of pre-inflationary dynamics and the comparison between the classical theory 
and LQC. 

It is known that due to the presence of anisotropic shear, the Hubble friction during cosmic 
expansion is enhanced which, in turn, leads to faster decay of the kinetic energy of the inflaton pQ]. 
This helps the isotropic slow-roll conditions arrive more quickly. If the inflaton is initially taken to 



roll down the potential at some initial time t = ti {(l){ti) < 0), then due to the enhanced Hubble 
friction, higher anisotropy leads to an increase in the number of e-foldings. In this article, in the 
discussion of classical dynamics, we consider a more generic initial condition on the initial velocity 
of the inflaton, including the case when the inflaton is initially rolling up the potential in the pre- 
inflationary era. We find that, for such an initial condition i.e. </>(ti) > 0, the number of e-foldings 
decrease with an increase in the anisotropic shear in the classical theory in comparison to the 
corresponding isotropic evolution. Further, in the regime of low anisotropy, there is a significant 
difference in the number of e-foldings depending on whether initial value of (j) is positive or negative. 
In contrast, we find that the large anisotropy gives rise to the same number of e-foldings for both 
the signs of initial (j). The latter behavior of the amount of inflation can be attributed to the 
strength of the Hubble friction which, at large anisotropy, is so strong that the sign of (p becomes 
insignificant. 

In the classical theory, for a given type of initial condition, the number of e-foldings (A^) show 
a monotonic variation with anisotropic shear i.e. A^ either decreases (for an inflaton rolling up the 
potential) or increases (for an inflaton rolling down the potential) with an increasing initial value 
of the anisotropic shear. In LQC, however, it turns out that the amount of inflation does not 
monotonically vary with the anisotropic shear. In the case when initial is chosen to be negative^ , 
the number of e-foldings increase with increasing shear if the anisotropic shear scalar cr^ is less than 
a particular value o"^*, and attain a maximum value when a^ = a'^^,. Thus, unlike the classical 
theory, there exists a range of a^: a'^^, < a'^ < cr^ax) fo^' which an increase in anisotropy decreases 
the number of e-foldings. On the other hand, if the initial velocity of the scalar field is positive, 
the behavior of N is opposite. That is, A^ decreases for a^ < cj^*, increases if o"^ > o"^* and attains 
a minimum value at o"^ = o"^*. Interestingly, the value of o"^* turns out be independent of whether 
the initial velocity of the scalar field is positive or negative (though it depends on the absolute 
values of initial (p and initial (p) . 

We also study the attractor behavior of the phase-space trajectories of the classical and LQC 
effective dynamics. In the classical theory, it is known that all the classical trajectories join the 
isotropic slow-roll inflationary trajectory in their future evolution provided the initial value of 
(f) is high enough to allow inflation (see for eg. [22]). Also, in the limit of initial shear scalar 
tending to zero, Bianchi-I trajectories approach the isotropic Friedmann-Robertson- Walker (FRW) 
trajectory, irrespective of the initial energy density. As in GR, it turns out that in LQC, irrespective 
of the initial anisotropic content of the spacetime, all the Bianchi-I spacetime trajectories, meet 
the isotropic slow-roll inflation in their future evolution, if the initial value of the inflaton field is 
high enough. In this way, slow-roll inflation is an attractor for all the trajectories in the effective 
spacetime of Bianchi-I LQC with a (j)'^ potential starting at the bounce. Unlike in the classical 
theory where there exists an isotropic trajectory for every given value of energy density in the 
initial data, in LQC, the energy density at the bounce is fixed at p = pmax in the isotropic 
spacetime. Therefore, in order to obtain the isotropic limit of Bianchi-I spacetime in LQC at the 
bounce, one has to consider p — )• /Omax in addition to cr^ — )• 0. If the energy density at the bounce 
is fixed at any value other than Pmax (which is allowed in Bianchi-I LQC), then it turns out, that 
the shear scalar can not be decreased to zero. In that case, there is a minimum non-zero value of 
the shear scalar, depending on the energy density at the bounce. In this sense, the approach to the 
isotropic limit in the effective description of LQC is subtle in comparison to the classical theory. 
That is, in the classical theory, isotropic spacetime can be approached by decreasing a'^ to zero for 
any fixed value of energy density in the initial data, whereas, in LQC, the energy density at the 



^ The initial conditions in most of the simulations for LQC discussed in this work are provided at the bounce of the 
mean scale factor. 



bounce must also tend to Pmax in addition to o"^ — t- 0. 

This manuscript is organized as follows, in section-II, we present the dynamical equations of 
Bianchi-I spacetime in the classical and the effective description of LQC in the Ashtekar variables. 
We discuss the features of anisotropic geometries, slow roll parameters, conditions for accelerated 
expansion, and behavior of equation of state of the scalar field in the classical theory. Based on 
the modified dynamical equations we discuss some key properties of the mean Hubble rate and the 
shear scalar in the effective description of LQC. In section-Ill, we discuss the numerical studies 
of the model where we present the results of various numerical simulations of classical and the 
LQC dynamics. In all the simulations, the mass of the inflaton is fixed to m = 1.21 x lO^^rripi, 



consistent with the WMAP data |65l |. With the help of explicit numerical simulations, we show 
the occurrence of non-singular bounces for various values of initial anisotropic shear in the effective 
dynamics of Bianchi-I spacetime. We explore the process of isotropization, behavior of the number 
of e-foldings for various initial conditions and study the phase space trajectories of the classical 
and LQC Bianchi-I spacetime. To discuss the qualitative behavior of trajectories, the initial value 
of inflaton in various plots is taken as (p = 3.14m,pi. Since this initial value is also required in 
the classical isotropic inflationary model for sufficient number of e-foldings, it serves as a good 
representative of the inflationary trajectories to study the effects of presence of anisotropy on the 
amount of inflation as compared to the classical isotropic slow-roll. Summary of results from 
various other initial values of the inflaton is presented in the tabular form, in Tables- HI and HIl We 
summarize with the discussion of results in Sec. V. 

II. CLASSICAL AND EFFECTIVE DYNAMICS: BASIC EQUATIONS 

In this section we present the description of Bianchi-I spacetime in terms of the loop quantum 
variables - the Ashtekar-Barbero connection A^ and the triads Ef. In the following, the first 
subsection is devoted to the classical theory. After discussing the relationship of connection and 
triad variables with the metric variables, we outline the derivation of classical dynamical equations 
from the Hamiltonian constraint, and discuss in detail the conditions for accelerated expansion in 
the classical Bianchi-I spacetime sourced with a m'^(p'^ inflationary potential (with m = 1.21 x 10^^ 
mpi). In the second subsection we outline the derivation of effective dynamical equations in the 
effective spacetime description of Bianchi-I spacetime in LQC, and note key non-trivial features 
in comparison to the classical theory which play an important role in distinctions of dynamical 
behavior of trajectories in GR and LQC. 

A. Classical theory 

We consider a homogeneous Bianchi-I spacetime with the spatial manifold having the topology 
S(= M'^), thus endowing S x M topology to the spacetime manifold. The metric for Bianchi-I 
spacetime is given as 

ds^ = -N^ df + al dx"^ + aj dy^ + al dz^, (2.1) 

where aj denote the directional scale factors. In order to define a symplectic structure on the 
manifold we need to introduce a fiducial cell V on the spatial manifold of the spacetime, with 
fiducial volume Vq = hhh- Here k denote the co-ordinate lengths in the three spatial directions. 
The edges of the fiducial cell are chosen to lie along the fiducial triads e^. The fiducial metric 
compatible with the fiducial co-triads oj^ is denoted as Qab- 

Utilizing the underlying symmetries of the homogeneous spatial manifold of Bianchi-I spacetime, 
the connection A^^ and triads Ef can be written in terms of the symmetry reduced connections (cj) 



and triads (pi) as 

Ai = c\hr'u;i and E^ = Pi{k)V-' y^, (2.2) 

where q = det (qab) is the determinant of the fiducial metric and the index i runs from 1 to 3. 
The symmetry reduced connection Cj and triad pj are conjugate to each other and they satisfy the 
following Poisson bracket relation 

{c\pj} = 87rGj5i, (2.3) 

where 7 = 0.2375 is the Barbero-Immirzi parameter whose value is fixed using the black-hole 
entropy computation in LQG. The triads pi are kinematically related to the usual metric variables 
of the spacetime in the following manner, 

Pi = £ihh\a2a3\, P2 = £2hh\aia3\, p3 = e3^i^2|aia2| (2.4) 

where e, = ±1 depending on whether the triads have positive or negative orientation. In the 
following, without any loss of generality, we choose e^ = +1, and the coordinate lengths li to be 
unity. 

1. Dynamical equations 

The classical Hamiltonian constraint for Bianchi-I spacetime, with lapse A^ = 1, in terms of the 
connections Cj and triads pj can be written as follows [40] 

^ci = - „ ry ■?jr (ciC2PiP2 + cyclic terms) + ?^matt, (2.5) 



where V = y/pvPiPi stands for the physical volume of the cell V and "Hmatt denotes the Hamiltonian 



of the matter sector. Considering the form of the self interacting potential given via l^((/>) = 2'm?(l>^, 



the matter Hamiltonian of the system can be written as 

^matt = „ , '^ + -mVVl'iPzPs, (2.6) 

^^JPlP2P'i 2 

where P(f, is the conjugate momentum of the scalar field (p. Starting from the total classical 
Hamiltonian constraint, the dynamical equations can be derived using Hamilton's equations of 
motion 

dUci 



Pi = {pi^Tici] = -SnG-f- 



da 



Ci = {q,?^el} = 8^G7^. (2.7) 

OPi 

Further, the directional Hubble rates are given in terms of the time derivatives of the triads pi as 
follows: 

2 \P2 P3 Pi 

and similarly for H2 and H3 . Using the equations of motion and the expression for the directional 
Hubble rates it is now straightforward to write the classical Einstein's equations for Bianchi-I 
spacetime as 

H1H2 + H2H3 + H3H1 = 47rG(</.2 + m202), 
H2 + H3 + HI + HI + HI = -47rG(</)2 - m2</.2), (2.9) 



and the cyclic permutation of these two equations. Equations of motion for the scalar field can 
also be derived in a similar fashion which yield the Klein-Gordon equation, 

4> + 5H<j) + V^ = , (2.10) 

where H = {Hi+H2 + H-i)/3 is the mean Hubble rate and V^tf, is the derivative of the self interacting 
potential with respect to the field (j). The Klein-Gordon equation governing the evolution of the 
scalar field is equivalent to the conservation equation p = —3H{p + P), where p = (/>^/2 + V{(l)) 
is the energy density and P = 0^/2 — V{(^) is the pressure of the scalar field. In the inflationary 
scenario, the term 3H(f), also known as Hubble friction, plays an important role by causing the 
slow-roll of the inflaton in an expanding universe. A bigger friction term results in faster decay of 
the kinetic energy of the inflaton thus leading to an early start of the slow-roll conditions. This 
turns out to be the key to understanding the effect of the anisotropy on the amount of inflation 
obtained, as the anisotropic shear interacts with the dynamics of the scalar field through eq. (j2.10p . 
The anisotropic shear is measured by a scalar o"^ which in terms of the directional Hubble rates 
is given as 

,2 _ f^u. zT_^2 I f u_ t7'_^2 I / TT_ u \2 



a' = - i^{Hi - H2Y + {H2 - H^Y + (^3 - ^i)'j ■ (2.11) 

Using the expression of shear scalar o"^, the vanishing of the Hamiltonian constraint T-Ld ~ 0, results 
in the generalized Priedmann equation for diagonal Bianchi-I spacetime: 

^^9 SvrG In , ^ 

H^ = —-p+-a^. (2.12) 

Further, using Einstein's equations for diagonal Bianchi-I spacetime (given by eq. (j2.9p ) one can 
derive the Raychaudhuri equation as follows 

«=-i^(p + 3P)-l.^. (2.13) 

ad 6 

Eq. (|2.13|) can also be obtained by computing the time derivative of the Friedmann equation given 
by eg. (j2.12p and utilizing the conservation equation p = —3H{p + P) along with a'^ = —6Ha^ 
[6 a . |67t] . Note that in the classical theory, a"^ oc a~^ . 

2. Kasner exponents and Jacobs ' parameter 

To understand the structure of spacetime in the very early stage of evolution (including the 
pre-inflationary epoch), it is convenient to use the Kasner exponents, ki, which are related to 
the directional scale factors via Oj oc t '. The Kasner exponents have been extensively used to 
understand structure of the spatial geometry during the approach to classical singularity, which 
can be classifled into point, barrel, pancake and cigar types. When all the Kasner exponents are 
positive {ki,k2,k^ > 0), then all the directional scale factors approach singularity together and 
the spatial geometry of the spacetime tends to a point like structure. Similarly, the approach to 
singularity is pancake when two of the Kasner exponents are zero and one of them is positive 
(fci, /c2 = and A;3 > 0), barrel, when one of the exponents is zero while other two positive {ki = 0, 
^2, ^3 > 0), and cigar, when one of the exponents is negative and two of them are positive {ki < 
and ^2 ) ^3 > 0) . 

Using the deflnition of the shear scalar (cr^) and the expansion scalar 9 = Hi + H2 + -^3, one 
can write the quantity a'^/9'^ in terms of the Kasner exponents as follows 

cr^ 2 

■^ = -{l-3{kik2 + k2k3 + kski)). (2.14) 



The ratio of the shear scalar and the expansion scalar takes the minimum value when all of the 
three directions of the spacetime are expanding i.e. all the three Kasner exponents are positive. It 
can be easily proved that if a^ jQ"^ < 1/6, then the only solution possible is that of all of the Kasner 
exponents are positive. This corresponds to a point-like approach to the classical singularity in 
the backward evolution. Similarly, pancake type structure is characterized by o^ jQ'^ = 2/3, the 
cigar type structure takes place when a^ jQ"^ > 1/6 and the barrel type structure is formed when 
1/6 < G^ jQ'^ < 2/3. Thus, a'^ /O"^ describes a hierarchy of the spatial geometry based on the 
anisotropic shear present in the spacetime [68|. Since the expression for a'^/6'^ does not assume 
any particular form of matter, the spatial structure of the geometry given by eq. (I2.14p holds true 
irrespective of the matter content. As an example, for the inflationary spacetime, if a'^/6'^ < 1/6 
in the future evolution, then all the Kasner exponent must take positive values and hence all the 
directions of the spacetime must be expanding. 

For the simplicity of the following calculations, we use the Jacobs' parameter ej = (T/(47rGp)^'^ 



671 ]. The parameter ej denotes the dominance of the anisotropic shear over the matter density. 
A large value of ej refers to a highly anisotropic universe while ej — >■ corresponds to isotropiza- 
tion (i.e. when the contribution of matter density to the spacetime curvature is far greater than 
that of the anisotropic shear). Moreover, the value of ej can be related to the structure of the 
spatial geometry near the classical singularity. As an example, for ej < 2/v3 (which corresponds 
to a'^/6'^ < 1/6), all the directional Hubble rates are necessarily positive and the approach to 
singularity in the backward evolution is point like. For very large anisotropy, with ej ^> 2/v3, the 
approach to singularity is cigar like. For a detailed discussion on the relation between anisotropy 
and the hierarchy of the approach to classical singularity in diagonal Bianchi-I spacetime, see Ref. 



3. Accelerated expansion 

The condition for an accelerated expansion in terms of the equation of state and the parameter 
ej can be obtained by putting the condition a /a > in the Raychaudhuri equation (eq. (I2.13P ) for 
Bianchi-I spacetime. It yields: 

w<--^, (2.15) 

where w = P/p is the varying equation of state of the scalar field. For ej = the above expression 
reduces to the condition for inflation in isotropic spacetime {w < —1/3). 

According to the inequality (I2.15p . whether or not a Bianchi-I universe undergoes acceler- 
ated expansion depends on the position of the field in the infiationary potential and also on the 
anisotropic shear in the universe. This is in contrast to the isotropic model, where the position 
of the infiaton in the infiationary potential is sufficient to determine whether the expansion of the 
universe accelerates or not. We demonstrate this in Fig.[Tl which shows two plots depicting the 
validity of the inequality given in eq. (j2.15p . In this figure, the initial values specified at t = are 
: i;^(0) = —2 X 10~^ rrip^ and (j){0) = 3.14 mpi. (The qualitative behavior of the plot remains the 
same for other values of initial (p and (p). We choose (/>(0) = 3.147Tipi, just as a representative value 
to discuss the behavior of inflationary trajectories, as this value is close to the value of the inflaton 
field near the onset of slow-roll in the isotropic spacetime, to generate approximately 60 e-foldings. 
In the rest of the article, unless a different initial condition is specified, we will take 0(0) = 3.14?TT,pi, 
wherever required to discuss qualitative behavior of amount of infiation and dynamical trajectories, 
shown in the plots in the following sections. The left plot in Fig.[T] is for a small value of initial 
anisotropic shear and the right one corresponds to a higher anisotropic shear. In these plots, the 
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FIG. 1; In this figure the (red) solid curve shows the equation of state, the (blue) dashed curve denotes 
— (1 + ej)/3, the thin black horizontal line shows the w = —1/3 line and the black dot marks the onset 
of accelerated expansion (inflation). The accelerated expansion starts when equation of state satisfies the 
condition given by eq. (|2.15p . The initial data for these plots are 0(0) — —2 x 10~^ rripj and (j){0) — 3.14 mpi, 
with ej(0) = 4.01 for the left, and ej(0) = 1.17 x 10^ for the right plot. The mass of the inflation is taken 
to be m = 1.21 x lO'^mpi. 



initial conditions are taken such that the kinetic energy is greater than potential energy and the 
equation of state is initially^ w ~ 1. As the further evolution takes place, the kinetic energy decays 
and the equation of state decreases. After some time, the potential energy dominates over the 
kinetic energy and the equation of state becomes w ~ —1. In the right plot of Fig.[Tl we see that 
though the equation of state w attains the value —1 quite early, around t ^ l,000tpi, but, the 
accelerated expansion does not start until eq. (j2.15p is satisfied, which occurs around i ~ 50, 000 ipi. 
In contrast, in the left plot since the initial anisotropic shear is weak, the accelerated expansion 
starts close to tt; < —1/3. 

Using eq. (|2.15p . one may be tempted to conclude that a non- vanishing value of the parameter ej 
is detrimental to the amount of inflation in comparison to the isotropic model, since more negative 
value of the equation of state is required for an accelerated expansion of the universe. However, 
as earlier pointed out in Ref. [20[, this expectation turns out to be incorrect. It turns out that 
the presence of anisotropy actually helps inflation by increasing the Hubble rate, which in turn 
enhances the friction term [2>H(j)) in the Klein-Gordon equation. Due to this enhanced Hubble 
friction, the kinetic energy of the scalar field decays at a faster rate. This brings the condition for 
accelerated expansion to be achieved earlier, in presence of anisotropy, causing a longer phase of 
accelerated expansion. 

The condition for accelerated expansion can also be written in terms of the contributions from 
kinetic and potential energies, and the shear scalar, by using eq. (j2.13p . Substituting the expressions 
for energy density and pressure in terms of the kinetic (</'^/2) and potential energy V{4)), eq. ()2.13p 
yields. 
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(2.16) 



^By starting with kinetic dominated initial conditions we make sure that we start witli conditions not favorable for 
inflation to occur immediately. The initial trajectory, although away from the inflationary trajectory, joins the 
inflationary trajectory in the forward evolution (as discussed later in this article). 
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FIG. 2: This figure shows the evolution of kinetic {(jP/2), potential (m^(/)^/2) and anisotropic shear 
(ct^/IGttG) terms with time. The (red) solid curve shows the kinetic energy contribution, the (blue) short 
dashed curve corresponds to potential energy contribution, and the (orange) dashed curve shows the shear 
term. The black dot marks the onset of inflation. Inflation takes place in the regime when potential energy 
dominates over the kinetic energy as well as the shear term. 



The condition for accelerated expansion can then be vi^ritten as: 



a 



''w>nt+i6.G 



This relation tells us that for inflation (accelerated expansion) to take place, the potential energy 
must dominate both the kinetic energy and the anisotropic shear term in the above equation. In 
Fig.[2l we show the evolution of the kinetic and potential energy, and cr^/lGvrG with respect to time. 
It is evident that inflation does not start (denoted by a black dot in Fig. [2]), until the sum of shear 
term cr^/lGvrG and the kinetic energy term becomes less than half of the potential term. We also 
see that following the onset of inflation, there is a short period where the shear term dominates 
over the kinetic energy, before it decays quickly during accelerated expansion. Thus, there may 
be a short duration of anisotropy during the classical inflationary phase of Bianchi-I spacetime. 
However, this is not a generic feature of evolution in Bianchi-I spacetime. If the initial conditions 
are such that the shear energy is less than the kinetic energy in the beginning, then there will not 
be any duration of such dominance of shear over kinetic energy. 

For the subsequent discussion of the numerical results, it is useful to define the slow-roll param- 
eters. In the classical isotropic inflationary spacetime, the slow-roll parameters e and t] are defined 
as 



3^ 



</>2 + 2V{(P) 



and T7 = 

Hfp 



(2.17) 



The slow-roll inflation is defined as the phase of evolution where e, |r/| ^ 1. Smallness of these 
parameters implies that the Hubble rate {H) varies very slowly. In the discussion of the numerical 
results, we consider above definition of slow-roll parameters in the inflationary Bianchi-I spacetimes. 
It is evident from the eq. (|2.17p that for small slow-roll parameters, e, |?7| <C 1, the value of 
both (j) and (p will be small. This implies that during the slow-roll phase, the inflaton rolls down 
the potential with a small and almost constant velocity. An example of the numerical evolution 
of e and 4>, in Bianchi-I spacetime, starting from (^(0) = 3.14mpi, (/)(0) = 0.00002 rripj and ej(0) = 
1.58 X 10^^, is shown in Fig.O We see that, once e <C 1, the slow-roll inflation sets in, and the field 
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FIG. 3: This figure shows an example of the evolution of the slow- roll parameter e and corresponding 
field velocity (f). The dashed black line in the left plot marks e « 0.008 and in the right plot it shows 
the corresponding value of the field velocity (f) k, 1.97 x lO^^TTipj. The initial conditions are taken as 
0(0) — 3.14 mpi, (/)(0) = 0.00002 TOpj and ej(0) — 1.58 x 10^"*. It is clearly seen that even for initial anisotropy 
as high as ej(0) = 1.58 x 10^*, the slow-roll conditions are met in the future evolution. 



velocity varies very slowly. For these initial conditions, during this phase, the value of velocity of 



the inflaton is 



1.97 X 10 



-7 



m 



PI 



and the slow-roll parameter is e ~ 0.008. 



B. LQC Effective equations 



Canonical quantization of cosmological models in LQC is based on techniques in LQC. The 
classical gravitational Hamiltonian constraint of a given model is first expressed in terms of the 
holonomies of connections ^4" and fluxes of the triads E'^, which are then expressed as appro- 
priate quantum operators. The quantum Hamiltonian constraint hence obtained gives rise to a 
non-singular quantum difference equation [26]. Based on the geometric formulation of quantum 
mechanics, an approximate continuous description of the underlying quantum geometry, known 
as the effective description of LQC, can be obtained through a faithful embedding of the classi- 
cal phase space into the quantum phase space of LQC [51(]. Through an appropriate choice of 
semiclassical states, those which correspond to a macroscopic universe at late times, an effective 
Hamiltonian constraint can be derived by calculatirigthe leading order corrections introduced by 



46( 1 ■ The effective Hamiltonian constraint 



the underlying discrete geometry of the spacetime [4f 
hence obtained gives rise to the modified equations of motion, which introduce corrections to the 
classical Friedmann and Raychaudhuri equations. Extensive numerical simulations for various mat- 
ter models, including the Bianchi-I anisotropic spacetimes, by considering the LQC evolution of a 
semi-classical states show that the effective theory agrees very well with the full quantum dynam- 
ics at almost all scaleSjincluding the deep Planck regime and the classical low curvature regime 

SS, S H S m, IS-H- 

The effective Hamiltonian constraint for the Bianchi-I spacetime, with lapse A^ = 1, is given via 



•^ The effective Hamiltonian considered here does not incorporate "inverse triad corrections" which are only mean- 
ingful in the models with a compact topology. As mentioned earlier, the topology for the Bianchi-I spacetime, 
assumed to be E x TZ, is non-compact in the present context. Even if one considers a compact topology, since the 
physical volume of the fiducial cells in the numerical simulations remain large compared to the Planck volume, the 
corrections due to the inverse triad terms would be neghgible. 
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^eff — — 



sin(/iici) sm /i2C2 \ , n, /o1e^ 
= = P1P2 + cyclic terms + fima.tt, (2-18) 



5— 1 z = yiP2 -r ^y^i^^ uciiiio I 



where V = ^JvvPWi is the physical volume, T^matt is the matter Hamiltonian taken to be same as 
given in eq. (j2.6|) and /xj are lengths of the plaquette around which the holonomies are calculated 

/ii = A./-^, iL2 = \J-^-^ and /is = A./^-^ . (2.19) 



Pi V P2 V V-i 

Here A^ = 4\/37r7Zp[ denotes the minimum area gap which results from the quantum discreteness 
of spatial geometry in LQG. 

The effective Hamiltonian constraint given via eq. (I2.18P generates the dynamics of the spacetime 
for a given matter part of the Hamiltonian constraint, T^matt- The equations of motion of the triads 
and the connections are given via the Hamilton's equation of motion for 'H.^i^ (see for eg. [50]): 

V\ - - - 

pi = — - (sin(/i3C3) + sin(/i2C2)) cos(/iiCi) (2.20) 

'/A 



and 

ci = 



— — [cos(/i2C2) (sin(/iici) + sin(/i3C3)) C2P2 + cos(/i3C3) (sin(/iici) + sin(/U2C2)) C3P3 
7AP1 

- COs(//iCi) (Sin(//2C2) + Sin(/X3C3)) cipi 

BT-I 

(sin(/i2C2) sin(/i3C3) + cyclic terms) 1 + 87rG7 — ""^ . (2-21) 

The equations of motion of other directional triads and connections can also be computed in a 
similar fashion. The equations of motion hence obtained, along with the initial data which satisfies 
the effective Hamiltonian constraint given by eq. ()2.18p . lead to a well posed initial value problem. 
The dynamical equations obtained from the effective Hamiltonian constraint modify the classical 
Einstein's field equations in a significant way. It is evident from eq. (|2.20p that Pi/pi never diverges 
during the evolution as the value of sin(/ijCj) andcos(/[tiCj) are bounded above by unity. This 
further leads to non-diverging curvature scalars of Bianchi-I spacetimes such as energy density {p) , 
shear (cj^) and expansion {9) scalars for the comoving observers which are bounded above by their 
respective maxima [53] • These values are given by: 

2.78 , 2 11-57 

Pmax ~ U.4ippi, ^max ~ —. , and CF^^^ K, — — . 

'pi tp[ 

Moreover, it has been shown that all the curvature invariants of Bianchi-I spacetime are bounded 
and all the strong singularities of Bianchi-I spacetime, for the matter satisfying null energy condi- 
tion, are resolved [56}]. The resolution of singularity results in a non-singular bounce in the evolution 
of Bianchi-I spacetime in the effective description of LQC. Far way from Planckian regime, where 
the spacetime curvature is small, the effective dynamical trajectory agrees well with the classical 
theory. In the backward evolution as the spacetime curvature grows, the effective trajectory differs 
from the classical trajectory. During the subsequent backward evolution the classical trajectory 



Sometimes in literature, one finds a different functional dependence of fii on triads. It turns out tfiat tfie form 
Hi{pi,P2,P3) given by tfie equation above is the only known choice where resulting physics is independent of the 
rescaling freedoms of the fiducial cell when the manifold is non-compact (see ref . [52| for details) . 
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goes into singularity, whereas the effective dynamical trajectory undergoes a non-singular bounce. 
In the effective dynamical evolution, the curvature scalars are always finite and bounded by their 
respective maxima given in the above equation. 

The finite and non-diverging values of the curvature scalars introduce interesting modifications 
in terms of the choice of initial conditions and the dynamics of the evolution of universe in the 
pre-inflationary regime in the Bianchi-I spacetime in LQC. There are two important modifications 
relevant to our discussion. 

• Unlike in the classical theory, it turns out that in LQC, the mean Hubble rate does not 
monotonically vary with the shear scalar. This behavior has been noted for the first time in 
the numerical simulations discussed in this paper^. Due to such a non-monotonic behavior 
of the Hubble rate as a function of the shear scalar, the Hubble friction in the Klein-Gordon 
equation is affected in a similar fashion. In other words, unlike in the classical theory, Hubble 
friction is not always enhanced by the anisotropic shear in LQC. 

• The shear scalar in LQC during the pre-inflationary regime behaves very differently than that 
in the classical theory. Unlike in the classical theory, where a'^ oc a~^, shear scalar does not 
have a monotonic dependence on the scale factor, near the bounce. There may be instances 
when, depending on the conditions at the bounce, the shear scalar starts out with a very 
small value, and then increases to reach a local maxima. On the subsequent evolution, the 
shear scalar would decrease and approach the classical value in the regime where quantum 
gravity corrections are weak. As we will see in next section, this has important implications 
in regards to the isotropic limit and the attractor behavior of Bianchi-I spacetime in the 
vicinity of the bounce. 

III. CLASSICAL AND EFFECTIVE DYNAMICS: NUMERICAL RESULTS 

In this section we present numerical studies of the Bianchi-I spacetime with an inflaton as a 
matter field with a self-interacting potential V{(/)) = m?<f)'^ /2. In the first subsection, we discuss the 
results of the classical theory and the following subsection is devoted to the numerical simulations 
of the effective spacetime description in LQC. We discuss important features of both the classical 
and the LQC descriptions of the Bianchi-I inflationary spacetime, such as the effect of anisotropic 
shear on the amount of inflation, the process of isotropization, the isotropic limit of Bianchi-I 
spacetime, and the attractor behavior. Over a hundred simulations were performed for both the 
classical theory and LQC by varying different parameters including the initial shear scalar, the 
initial value of the inflaton and the initial velocity of the inflaton. In all of the simulations, the 
mass of the inflaton is taken as m = 1.21 x 10~^ m-p\ and the numerical accuracy of the vanishing 
of the Hamiltonian constraint is of the order of 10^^'^. Note that there is a symmetry associated to 
the signs of the value of inflaton field and its time derivative. That is, the results remain invariant 
if the sign of both (p and (j) are changed simultaneously. In the numerical simulations performed 
in this article, we take the initial value of the inflaton field to be positive. For an infiaton which 
is initially rolling down, we take (/>(0) < and for an inflaton which is initially rolling up, we take 
(/>(0) > 0. 



* Due to unavailability of the modified generalized Friedmann equation for the Bianchi-I spacetime in LQC, it is 
difficult to obtain this result analytically. However, in the regime where anisotropic shear is weak, the modified 
Friedmann equation for Bianchi-I spacetime in LQC has been obtained and it contains o"* corrections [69|. These 
corrections can cause the Hubble rate to decrease if the anisotropic shear increases beyond a certain value. 
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A. Classical theory 

As discussed in the previous section, the dynamics of the classical Bianchi-I spacetime is gov- 
erned by the generalized Friedmann and Raychaudhuri equations (eq. (|2.12p and eq. (|2.13p ) from 
where it follows that specifying the energy density, pressure and the shear scalar at a given time 
provides a complete initial data. For a scalar field the pressure and energy density can be ex- 
pressed in terms of the value of the scalar field (j) and its time derivative (p. Thus, giving the 
initial value of the scalar field, its velocity and the shear scalar forms a complete set of informa- 
tion required to compute the state of the universe at any later time. We evolve the initial data 
provided in a pre- inflationary phase and study various properties of the dynamical trajectories 
including isotropization, amount of inflation and the phase-space portraits. As discussed in the 
previous section, inflation takes place when the potential energy dominates over the sum of kinetic 
energy term [(jp' /2) and the shear contribution ((j^/lGvrG). In the numerical simulations presented 
in this section we will, however, start with a kinetic and shear dominating initial conditions. This 
gives insights on the way conditions in favor of inflation are achieved, starting from comparatively 
unfavorable initial conditions. Moreover, it brings out important features about the behavior of 
the pre-inflationary dynamics and its dependence on various parameters. 

1. Isotropization: behavior of the directional scale factors 

We first discuss the process of isotropization of Bianchi-I spacetime, that is the way a Bianchi-I 
universe starting with highly anisotropic initial conditions at the bounce with at least one scale 
factor contracting turns into one with all-expanding-directions. A Bianchi-I universe can evolve 
from highly anisotropic initial conditions which, depending on the strength of anisotropy, entail to 
various possible geometrical structures of the spatial geometry such as a cigar, pancake etc, close to 
the classical singularity. The key question is the way an inflationary spacetime with all expanding 
directions emerges from such a state. 

Recall that the structure of the spatial geometry of the spacetime can be determined via the 
value of the Jacobs' parameter ej. For example, ej < 2/v3 is a sufficient condition for all the 
directions to expand, whereas, for ej > 2/v3, not all of the directions will be in expanding state 
(as shown in the left plot of Fig. |4] where one of the directions is contracting while the other two 
expand in the forward evolution). Thus, at the onset of inflation, depending on the value of ej, 
all the directions may either be expanding or not. From the equations of motion of Bianchi-I 
spacetime, it is straightforward to see that the anisotropic shear (cr^) in an expanding classical 
Bianchi-I spacetime always decreases, as o"^ oc a^^ . During inflation, the mean scale factor grows 
exponentially which leads to a rapid decrease in a^, and hence in the value of ej. Due to this 
reason, ej quickly becomes less than 2/\/3 even if it starts from a larger value. Thus, even if one 
starts with a large ej in the Bianchi-I spacetime, ej becomes less than 2/v3 either before, or shortly 
after, the start of inflation, and all the directional scale factors turn out to be in the expanding 
state. 

Fig. [3] shows the evolution of the three directional scale factors starting from a cigar like struc- 
ture, such that two of the scale factors are increasing while the remaining one decreases, as the 
universe expands. It is evident from the left plot of Fig.H] that in the forward evolution, the de- 
creasing scale factor undergoes a turn around shortly after the time of onset of inflation, t = Tqn 
(defined as the time where a/ a becomes positive). After the turn around, all the scale factors 
undergo expansion. The adjoining right plot shows the evolution of individual Kasner exponents. 
In the beginning, one of the Kasner exponents, corresponding to the contracting directional scale 
factor, is negative which becomes positive soon after the onset of inflation. 
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FIG. 4: The left plot shows the evolution of directional scale factors in classical Bianchi-I spacetime in the 
inflationary scenario. The initial data is such that backward evolution gives a cigar like singularity. It is 
evident that, in the forward evolution, the contracting direction undergoes a turn around near the onset 
of inflation and during most of the inflation all the scale factors expand. The adjoining plot shows the 
evolution of the corresponding Kasner exponents. Note that although u; = — 1 is attained earlier, inflation 
does not begin until t = Tqn due to the presence of anisotropy. (Ton corresponds to the time at the onset 
of inflation) . 

In a similar way, one can study the isotropization of the Bianchi-I spacetime from the barrel 
and pancake like structures. In these cases too, irrespective of the their initial behavior, all the 
directional scale factors enter expansion phase either before or soon after the onset of inflation. 
In this way, an inflationary Bianchi-I spacetime isotropizes, no matter what the pre-inflationary 
conditions are. Such a turn around of the scale factor is a feature of a generic initial data, except 
for those with point like structures for which all the three directions are already expanding. Similar 
turn around of directional scale factors and the evolution of Kasner exponents was also noted in 
Ref. [22]. 

2. Number of e-foldings and initial conditions 

The number of e-foldings measure the amount of inflation during the inflationary era of the 
universe. This number is defined as A'^ = ln(oj/aj), where aj corresponds to the value of the 
mean scale factor where accelerated expansion ends, and Oj denotes the value of the mean scale 
factor where the accelerated expansion begins. We performed many simulations to compute the 
number of e-foldings for various initial values of the anisotropic shear. Fig. [5] shows the variation 
of the number of e-foldings as a function of the initial anisotropy (measured by ej), for a scalar 
field having a negative initial velocity (i;^(0) < 0, corresponding to a rolling down inflaton). To 
discuss the qualitative behavior of the amount of inflation, in the plots shown in this article, we 
choose the initial value of the inflaton to be 0(0) = 3.14 mpi. As mentioned earlier, this is just 
a choice representing inflationary trajectories and the way anisotropic shear affects the sufficient 
number of e-foldings, as this value is close to the value of inflaton at the onset of isotropic slow- 
roll to generate approximately 60 e-foldings in the isotropic spacetime®. The qualitative behavior 
discussed in this section also holds true for other values of (p{0) for which inflation takes place. 
We will consider other values of (^(0) while discussing the quantitative details of the amount of 



^ In order to obtain approximately 60 e-foldings in the classical isotropic spacetime, the slow-roll has to start at 
« 3.14 771P1 with the required value of the velocity of the inflaton being cf> « 1.97 x 10~ 



-7 2 

mpi 



14 



e-foldings. In confirmation with the known results in classical theory in this scenario, we find that 
if the inflaton is initially rolling down then the number of e-foldings increases with an increase in 
anisotropy. Further, different initial (/>'s produce the similar number of e-foldings for large initial 
anisotropy. That is, the spacetime tends to achieve the same number of e-foldings for large value 
of initial anisotropic shear irrespective of the initial kinetic energy of the scalar field. Note that for 
the higher values of initial anisotropy, an inflaton which rolls down the potential without initially 
satisfying slow roll conditions, quickly enters into slow roll. The universe starts inflating almost 
exponentially, which shortens the stage of anisotropic evolution. Numerical simulations show that 
for high values of initial anisotropic shear, the phase of anisotropic inflation is so short that any 
change in anisotropy produces little effect on the number of e-foldings. 
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FIG. 5: The left plot in this figure shows the number of e-foldings, N plotted against the initial value 
of In(ej) for various values of the initial 0. For each initial value of 0, the energy density is kept fixed 
while the anisotropic shear is varied. From top to bottom, curves correspond to 0(0) = (—2 x 10^^, —2 x 
10"'^, —2 X 10~^, — 0.1) TTipj. For large values of initial ej, all the trajectories tend to produce the same 
number of e-foldings. The right plot shows the variation of the number of e-foldings with time for various 
initial anisotropy parameters ej, but starting with the same initial value of (p: (f>{0) = —0.1 mpi. In both the 
plots, the initial value of the inflaton field is taken as 0(0) = 3.14 mpi. 



The right plot in Fig. [5] shows the evolution of the number of e-foldings during the forward 
evolution of inflationary Bianchi-I spacetime starting with the same 0(0) = — O.lmpj but for the 
various values of the initial anisotropy parameters ej. For these initial conditions, the number of e- 
foldings in the isotropic case (the case ej = 0) is less than 10. Whereas, in the Bianchi-I spacetime, 
the same initial conditions on the scalar field produce significantly more number of e-foldings for 
higher ej. Thus for 0(0) < 0, anisotropic shear not only increases the number of e-foldings, but a 
large anisotropy also makes slow-roll conditions be achieved irrespective of the initial velocity of 
the inflaton due to an increase in the Hubble friction in the Klein-Gordon equation.^ 

We now consider the case of inflaton initially rolling up the potential, i.e. 0(0) > 0. To our 
knowledge, this case has not been discussed elsewhere in the literature. The left plot in the Fig. [6] 
shows the variation of number of e-foldings with the increasing anisotropy in the initial data for 
various values of 0(0) = (2 x 10~^, 2 x 10~^, 2 x 10~^, 0.1) mpj. In all these simulations, the value 
of the inflaton field in the initial data has been taken as 0(0) = 3.14?7ipi, on the same footing as 
in Fig.O In contrast to the case of inflaton rolling down the potential, it turns out that, in this 
case, the number of e-foldings decrease as the anisotropy in the initial data increases. For large 



In GR, irrespective of the initial inflaton velocity, one can always find a suitable value of initial anisotropic shear 
such that slow-roll is obtained. 
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FIG. 6: This figure shows the number of e-foldings, N plotted against In(ej) for various values of the initial 
^. In the left plot, from bottom to top, the curves correspond to 0(0) = (2 x 10~^, 2 x 10"'^, 2 x 10~^, 0.1) mpj 
and initial value of field is taken 0(0) = 3.14 mpi for all the simulations. The right plot shows the time 
evolution of the number of efoldings for various initial anisotropy starting with the same initial value of the 
inflaton velocity 0(0) = 0.1 mp\. 



initial anisotropic shear, all the trajectories attain the same number of e-foldings irrespective of 
the initial (/). 

To understand the decrease in number of e-foldings with an increase in anisotropy, let us for a 
moment consider the isotropic model. In the isotropic case, an initial 0(0) > causes the field to 
roll up the potential before the onset of inflation. For example, if the initial 0> is say 0(0) = 0.1 m,pj 
(with initial as 3.14mpi), then the inflaton will roll up the potential and inflation starts to take 
place at ~ 4.8 mpi. As a result, the amount of inflation is greater as compared to the case 
of initial < 0. On the other hand, in Bianchi-I spacetime due the increase in Hubble friction, 
there is a faster decay of the kinetic energy of the field, and the field stops its upward journey on 
the potential before reaching ~ 4.8 mpi. Thus, infiation starts taking place at a comparatively 
smaller value of the field. Hence, for an inflaton which is initially rolling up the potential, an 
increase in anisotropy reduces the the number of e-foldings in comparison to the corresponding 
isotropic evolution. At very large anisotropy, the Hubble friction is so strong that the kinetic 
energy of the field decays quickly, and inflation onsets very close to the initial value of the inflaton. 

It is important to note that for a large initial anisotropy, both the positive and negative initial 
velocity of the field give the same amount of infiation. This is evident from the comparison of Fig. [5] 
with Fig.m It is also clear from these plots, that although large anisotropic shear may diminish the 
number of e-foldings, yet, it never forbids the occurrence of inflation. The right plot in Fig.[6l shows 
the time evolutions of the number of e-foldings, for 0(0) > 0, starting with different anisotropy 
but with the same initial 0(0) = O.lrrapj. It is evident that the maximum number of e-foldings 
decreases with an increasing initial anisotropy. In Table-|T]we provide a representative summary for 
various simulations to determine the number of e-foldings obtained during the inflationary phase 
for various initial values of the 0(0) = (—0.002, 0.002, —0.02, 0.02)?7ip[, and different values of 
initial anisotropic shear for 0(0) = (1.0, 2.0, 3.0, 3.14, 4.0, 5.0) mpi. It is clearly seen that at large 
anisotropic shear, the number of e-foldings for a given 0(0) is almost the same for all 0(0). 



3. Phase portrait and attractor behavior 

We now discuss the phase-space trajectories of an infiationary Bianchi-I spacetime in the classi- 
cal theory. Fig. [7] shows the evolution of various classical trajectories in the 0-0 space for Bianchi-I 
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TABLE I: This table summarizes the number of e-folding in the classical Bianchi-I spacetime for various 
initial values of (j){0) and (j){Q), depending on value of anisotropic shear present in the initial data. Note that 
for Ej = 0, we have not shown the entries corresponding to (f> — 3.14r7ipi leading to 60 e-foldings. For the 
conditions required to obtain 60 e-foldings in the isotropic case, see footnote-E] in Sec.-IIIA.2 

spacetime. Trajectories in the same quadrant correspond to different initial anisotropics but with 
the same initial values of (p and (p, that is, the initial energy density is fixed for all the trajec- 
tories and initial shear scalar is varied. The horizontal (red) arrows show the order of decreasing 
anisotropy. The trajectories with labels and (2) correspond to the isotropic spacetime (i.e. 
ej = 0). Each isotropic curve is composed of two distinct phases. The vertical part (labeled by (T)) 
corresponds to the phase in which there is no slow-roll. The horizontal part (labeled by @) corre- 
sponds to the slow-roll phase in the isotropic case. The black dots on the trajectories denote the 
point when the accelerated expansion begins. The trajectories in the upper half of the phase space 
correspond to 4>(0) > and the curves in the lower half show the trajectories for 0(0) < 0. For all 
the curves, shown in the Fig. [71 initial conditions on the inflaton are kept fixed at \4>\ = 3.14mpi 
and |(;/'(0)| = 0.1 m|[, and the anisotropic shear is varied. It is clearly seen from Fig. [71 that all the 
trajectories corresponding to various initial values of anisotropic shear meet the slow-roll curve in 
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FIG. 7: This figure sliows tlie two attractor behavior of the infiationary Bianchi-I spacetime. The dashed 
trajectories (denoted by 0) show the isotropic non-siow roli FRW attractor while the horizontal lines to 
which all the trajectories meet depict the isotropic slow-roll attractor (denoted by @)- The directions of 
(red) horizontal arrows show decreasing anisotropy (i.e. increasing parameter ej). The black dots denote 
the onset of accelerated expansion, k is k — SttG/S. 

their evolution. In this sense, the isotropic slow-roll is an attractor for all such solutions. 

For (j){0), (j){0) > (also for (j)(0) , (j){0) < 0), i.e. the initial conditions corresponding to the 
inflaton rolling up the potential, trajectories in Bianchi-I model remain close to isotropic slow 
roll curve for a longer time for smaller initial value of the anisotropic shear. Furthermore, with 
increasing anisotropic shear, the duration of the slow-roll decreases. This behavior is compatible 
with the variation in the amount of inflation with the varying initial anisotropy discussed earlier. 
On the other hand, trajectories corresponding to (j){0) > and (j){0) < (or with (p{0) < and 
0(0) > 0) tend to have longer duration of slow-roll for larger value of the anisotropy, which is 
again compatible with the observation that an increase in anisotropy increases the e-foldings for 
an inflaton initially rolling down the potential. 

Fig. [8] shows the 3D plot of the phase space trajectories of the field in the classical Bianchi-I 
spacetime. It turns out that as the shear scalar tends to decrease, the phase-space trajectories in 
classical Bianchi-I spacetime approach that of the corresponding isotropic classical FRW spacetime 
(with o"^ = 0). This behavior is also evident from the Fig. [7] where the thick (red) arrows show 
the order of decreasing anisotropy. In this sense, in addition to the attraction behavior towards 
the slow roll trajectory, the non slow roll isotropic curve also acts as an attractor for trajectories 
in Bianchi-I spacetime in the limit o"^ — )• 0.^ In other words, the classical FRW spacetime is 
an isotropic attractor of the classical Bianchi-I spacetime. That is, as the initial anisotropy is 
decreased, the trajectories tend closer to the classical FRW isotropic trajectories (shown by the 



For an earlier discussion of the second attractor, see Ref. ^ 
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FIG. 8: This figure shows the 3D plot of the phase-space trajectories of the scalar field. The trajectories 
starting in the same quadrant of </) — plane correspond to the same initial values of and (j) but different 
initial shear scalar. The dashed trajectories in the cr^ = plane depicts the corresponding isotropic trajec- 
tories. It is clear that as the shear scalar is decreased the Bianchi-I trajectories tend to the corresponding 
isotropic trajectory. 

dashed-thick lines in Fig. [7] and Fig. [8]). To differentiate this attractor from the slow-roll, we call 
it the 'non slow-roll FRW attractor'. In this precise sense, there is a double attractor behavior in 
the classical Bianchi-I spacetime. Before we go into the discussion of numerical results for LQC, a 
remark about the two attractors is in order. 

Remark: The two attractors discussed above, namely, the non slow-roll FRW attractor and the 
isotropic slow-roll attractor, are two distinct properties of the Bianchi-I spacetime. The approach 
to the isotropic slow-roll is a property of the Bianchi-I spacetime with an inflationary potential. It 
is achieved in the future evolution when the slow-roll conditions are met, whereas, the non slow-roll 
FRW attractor can be approached in the pre slow-roll phase of the evolution when the value of the 
initial anisotropic shear is decreased. In contrast, the slow-roll attractor is approached only after 
the onset of inflation. 



B. Effective loop quantum dynamics 

In this subsection we discuss the results from numerical simulations of the effective dynamical 
equations of Bianchi-I spacetime in LQC with a m'^c/}'^ potential. The mass of the inflaton is taken 
to be same as in previous sections: m = 1.21 x W~^ mpi. Unlike the classical theory, LQC is 
devoid of initial singularity. The mean scale factor instead of vanishing in the backward evolution, 
undergoes a non-singular bounce. Hence, LQC provides a way to investigate the inflationary 
scenario in Bianchi-I spacetime in a non-singular setting. In this article we present the first such 
study of inflationary Bianchi-I spacetime in a non-singular scenario. The non-singular evolution, in 
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the effective description of LQC allows us to extend the study of dynamics of Bianchi-I spacetime 
prior to the onset of inflation, in the deep Planck regime. 

As discussed in the previous section, irrespective of the initial conditions, all the curvature 
invariants of Bianchi-I spacetime are bounded and never diverge during the evolution given by the 
effective dynamical equations. This implies that the occurrence of bounce is a generic feature of the 
effective description of Bianchi-I spacetime in LQC, and bounce is obtained without violating any 
energy condition. If one starts with a classical expanding universe, so that the spacetime curvature 
is too small for the quantum geometric effects to be important, then during the backward evolution, 
initially the classical trajectory agrees well with the effective trajectory. In the further backward 
evolution, as the curvature of spacetime grows, the quantum geometric effects start becoming more 
prominent. Due to this, the effective trajectory starts deviating from the classical trajectory. In 
the subsequent backward evolution, the classical trajectory meets the initial singularity, when the 
mean scale factor goes to zero. On the other hand, the effective trajectory undergoes a smooth 
non-singular evolution, with the mean scale factor bouncing from a finite non-vanishing value. 
The individual scale factors, also, undergo smooth evolution from one side of the bounce of the 
mean scale factor to the other side. Figs. dH^a) and [TOlfa)) show two examples of evolution of 
the mean and the directional scale factors of Bianchi-I spacetime in LQC, with initial conditions 
being provided far from the bounce, close to the onset of inflation. The corresponding evolutions 
of the energy density and shear scalar are shown in Figs. ^Hh, c) and [TOlfb. c)). The initial 
conditions for the evolutions shown in Fig. [9] are 0(0) = 3.14?7ipi, 0(0) = —1.97 x 10^"^ rrip^, 
o"^(0) = 1.09x 10^^^ Ip^, and the initial mean scale factor is a(0) = 8.02x 10'^ while the initial energy 
density being p{0) = 7.23 x 10~^^ rrapj. Similarly the initial conditions for the evolution shown in 
Fig.Hniare 0(0) = 3.14 mpi, 0(0) = -1.99 x 10-^m|i, cj2(0) = 9.32 x W^^lpy, a(0) = 8.05 x 10^ 
and p(0) = 7.24 x lO^^^rripj. It evident from Fig. ([9l and [TO]) that the energy density and the 
shear scalars do not diverge, and are below their maximum values Pmax and cr^ax throughout the 
non-singular evolution. The values of the energy density and shear scalar at the bounce can be 
quite different, as can be seen by comparing Figl9] and FigllOi It is worth emphasizing that the 
bounce occurs without any fine tuning of the initial conditions or parameters. 

Since the bounce is a generic property of all the solutions in the effective dynamics of LQC, it 
provides a natural point in the evolution to specify initial data. In all the simulations which we 
discuss in the following, the initial conditions are provided at the bounce, characterized by i = 0. 
The bounce is characterized by the turn around of the mean scale factor which entails vanishing 
mean Hubble rate, H{0) = at the bounce. Therefore, the initial conditions along with satisfying 
the Hamiltonian constraint, must be chosen such that the mean Hubble rate equals zero. To 
ensure this, we provide pi{0), P2(0), P3(0), ci(0), 0(0), and 0(0) at the bounce, and compute the 
value of C2(0) and 03(0) by requiring that both the Hamiltonian constraint and the mean Hubble 
rate vanish. Keeping Pi{0), 0(0), 0(0) fixed, and by varying ci(0) at the bounce, we can vary the 
initial anisotropic shear scalar at the bounce for a fixed value of the energy density. This gives 
rise to a continuous set of solutions with a fixed energy density and varying initial shear scalar, all 
starting at the bounce. These solutions turn out to be the key to understand the effect of varying 
anisotropic shear scalar on the number of e-foldings. To illustrate. Fig. [11] shows the bounce of 
mean scale factors for varying (7^(0) for the same 0(0) = 3.14rn-pi and 0(0) = — 0.002 m,pj at the 
bounce. As in the previous sections, we choose the initial value of the inflaton field same as above to 
discuss the qualitative behaviors shown in the plots in this subsection. Results for different initial 
values of inflaton fleld will be given in the tabular form (see Table-lH|) . Another way of generating 
initial conditions with varying initial shear scalar is to vary the initial energy density by varying 
the initial velocity of the inflaton. This method of choosing initial data is essential in order to 
understand the isotropic attractor behavior of Bianchi-I spacetime in LQC, as we will see later in 
this subsection. 
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FIG. 9: The figure (a) shows the evolution of the mean and directional scale factors, close to the bounce, 
when initial conditions are given far from the bounce near the onset of inflation. Figures (b) and (c) show the 
evolution of energy density and the shear scalar respectively. The values of p and a^ are in Planck units. The 
initial conditions for this figure are 0(0) w 3.14 mpi, 0(0) w -1.97 x lO-'^mlj and cr2(o) ^ 1.09 x lO^^^^pi^. 
It is clear from the plots that the mean scale factor bounces while the directional scale factors undergo 
smooth evolution. 

In the following, we first study the evolution of the individual directional scale factors and the 
isotropization of Bianchi-I spacetime as the future evolution takes place. Then, we turn to compute 
the amount of inflation in terms of the number of e-foldings generated during the inflationary era, 
for a variety of initial conditions including different initial (j), (j) and cr^ at the bounce. We also 
study the evolution of dynamical trajectories in the <J3 — 4> phase space, starting from different initial 
conditions. The study of the phase trajectories brings out important features of the attractor 
behavior of Bianchi-I spacetime in effective description of LQC. 

1. Directional scale factors and isotropization 

We study the evolution of the spacetime starting from initial conditions characterized by the 
presence of a significant anisotropic shear. Due to the non- vanishing anisotropic shear, as in the 
classical theory, depending on the strength of the anisotropy, the expansion of the spacetime can be 
highly anisotropic. That is, while the mean scale factor expands, one or two of the directional scale 
factors decrease and the remaining one(s) increase in the forward evolution. However, in contrast 
to the classical theory, the backward evolution in LQC is non-singular. Here, we are interested in 
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FIG. 10: The figure (a) shows the evolution of the mean and directional scale factors, close to the bounce, 
when initial conditions are given far from the bounce near the onset of inflation. Figures (b) and (c) show the 
evolution of energy density and the shear scalar respectively. The values of p and a^ are in Planck units. The 
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initial conditions for this figure are (/)(0) sa 3.14topi, 0(0) « -1.99 x lO^^m^j and cr^(O) = 9.32 x 10"^"; 
It is clear from the plots that the mean scale factor bounces while the directional scale factors undergo 
smooth evolution. 



investigating the way an anisotropically expanding universe at the bounce turns into isotropically 
inflating spacetime in LQC. 

In Fig. 1121 we show the evolution of scale factors in a representative simulation in the Bianchi- 
I spacetime in LQC. In the expanding branch of the mean scale factor a, two scale factors are 
expanding and one is contracting near the bounce. This corresponds to a cigar like structure of the 
spatial geometry close to the bounce, in the backward evolution of the expanding branch. It turns 
out that long before the onset of accelerated expansion, LQC trajectory enters into the classical 
domain, and the shear scalar monotonically falls as ~ a~^ . Due to the ever decreasing shear scalar, 
in this domain, the parameter ej keeps decreasing and, just like in the classical theory, reaches 
ej < 2/\/3 before or right after the onset of inflation. As the accelerated expansion starts to take 
place, like in the classical theory, the parameter ej falls even more quickly and the contracting 
directional scale factor turns around, causing isotropization of the Bianchi-I universe in LQC. 

In Sec. IIA3, we discussed the evolution of kinetic and potential energies and compared them 
with the anisotropic shear contribution (u^/IG-ttG) in the classical theory. We now study their 
behavior in LQC. Fig. 1131 shows an example of the evolution of all of the three in an expanding 
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FIG. 11: This figure shows the evolution of the mean scale factor for different initial shear scalars (but with 
the same initial (p — 3.147tt,pi and (p — — 0.002 rripj) given at the bounce. The values of the shear scalar 
shown in the figure correspond to the values at the bounce (in units of lp{)- 





5x10^ 1x10^ 



FIG. 12: This figure shows the evolution of directional scale factors for cr^(O) — 3A7 Ip^ at the bounce with 
initial (p — 3.14 mpi and (p = — 0.002 rripj. The right plot shows the isotropization of the directional scale 
factors. At the bounce one of the scale factors start out in decreasing state but as the inflationary era 
progresses, all the three directions turn to expand. 

branch of inflationary Bianchi-I spacetime. Starting with a shear dominant initial conditions such 
that kinetic energy dominates over the potential energy, the kinetic energy monotonically decays, 
whereas, the shear energy shows a non-monotonic behavior in the vicinity of the bounce i.e. it 
attains a local maximum and then continues to decrease. In the forward evolution, as the spacetime 
curvature approaches the classical conditions, the anisotropic shear decreases monotonically. This 
behavior of the anisotropic shear stands in contrast with the classical evolution where the shear 
term a^ monotonically decreases in an expanding universe for generic initial conditions. This 
brings out a key difference between the process of isotropization in classical and LQC evolutions. 
It is also clear from the plot that the onset of accelerated expansion (marked by the black 
dot) takes place when the potential energy wins over the shear and kinetic energy, just like in the 
classical theory. This is not surprising because, in the above numerical simulation, at the onset 
of inflation the spacetime is well approximated with a classical solution. Following the start of 
inflation there may be a period where the contribution of anisotropic shear is still greater than 
kinetic energy (while contribution of potential energy being larger than the both). However, the 
shear scalar quickly fades away giving rise to isotropization during the inflationary phase. In this 
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FIG. 13: Evolution of the kinetic {(jP /2), potential (m^0^/2) and anisotropic shear contribution (ct^/IGttG) 
in LQC effective dynamics. The short dashed (blue) curve represents potential energy, the dashed (orange) 
line corresponds to the shear contribution and the solid (red) line shows the evolution of kinetic energy of 
the inflaton. 

way, similar to the classical theory, there may be a short period in the inflationary era when the 
anisotropic shear dominates the kinetic energy, while the potential energy remains greater than 
both. 



2. Amount of inflation 



We discussed earlier that the amount of inflation in the classical theory increases with increasing 
initial value of the shear scalar, for an inflaton which is rolling down. Since the shear scalar increases 
the mean Hubble rate it enhances the damping term in the Klein-Gordon equation. This in turn 
causes the slow roll conditions to be achieved earlier. Unlike the classical theory, in LQC the 
mean Hubble rate is not necessarily a monotonic function of the shear scalar. Due to the higher 
order corrections introduced by the underlying quantum geometry, the modified form of generalized 
Friedmann equation for Bianchi-I spacetime may contain nonlinear terms^ in o"^. Fig. 1141 shows the 
variation of the mean Hubble rate with an increasing initial value of shear scalar. It is evident that 
the mean Hubble rate shows a non-monotonic behavior with increasing shear scalar, and attains 
a maximum at some value of o"^. A result of this is that, in LQC, in the presence of anisotropic 
shear the number of e-foldings can show a non-monotonic behavior with an increasing shear scalar 
(while keeping the (/>(0) fixed). Fig. [15] shows the variation of the number of e-foldings for various 
different values of the initial energy density and (/>(0) < at the bounce. It turns out that for 
small shear scalar, the number of e-foldings (A^) increases with the values of shear scalar, then N 
attains a maximum value at some a^ = a1. For any higher value of anisotropic shear than aj, 
the amount of e-foldings decreases. Note that, for the values of the parameters considered here, 
the changes in the number of e-foldings is very small, as can be seen in Fig. [151 In the numerical 
simulations, al turns out to be al ^ lA7lp^ for \4>{0)\ < 0.02 m|,j, and for |(^(0)| > 0.02 m^j it 
weakly depends on the initial data. At this point the occurrence of a^ is only a numerical result. 
Due to the unavailability of the generalized Friedmann equation for Bianchi-I spacetime in LQC, 
we do not yet have analytical argument to explain its existence. 



® For example, an approximate expression in the small shear limit, discussed in the Ref. 691], contains a^ terms 
with opposite sign to a^ term, suggesting non-monotonic behavior of the mean Hubble rate with shear scalar. 
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As in the classical theory, for an opposite sign of (j) in the initial data (in this case the initial data 
given at the bounce), the amount of inflation shows opposite behavior with the initial anisotropic 
shear. The numerical simulations performed for 0(0) > are shown in Fig.ll6i It is evident that 
for 0(0) > 0, in contrast to 0(0) < 0, A^ decreases with increasing o"^ in the small shear regime, 
while it increases for a"^ > al- Interestingly, the value of a1 turns out to be numerically same as 
that for 0(0) < 0. 

Such a variation of A^ with initial anisotropy is tied to the behavior of the mean Hubble rate. As 
discussed in the classical theory, an enhancement in the mean Hubble rate, results in an increment 
in the amount of inflation if the inflaton is initially rolling down, and a decrement if it is rolling up. 
In LQC, since the mean Hubble rate first increases, attains a maximum and then decreases with 
increasing initial shear scalar (Fig. ll4p . the amount of inflation varies in the same way for 0(0) < 
(Fig. lisp. That is, it first increases, attains a maximum and then falls down. The behavior of N 
for 0(0) > is opposite to that of 0(0) < 0. 
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FIG. 14: This figure shows the non-monotonic behavior of the mean Hubble rate at the onset of inflation 
with varying initial value of the shear scalar at the bounce. The left plot corresponds to a negative initial 
(j) and the right plot corresponds to a positive initial cf) (in Planck units). 
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FIG. 15: This figure shows the variation of the number of e- foldings with the changing shear scalar at the 
bounce for < at the bounce. The values of is given in the units of rn^^. 

So far we have discussed the qualitative behavior of the amount of inflation with varying initial 
anisotropy at the bounce, let us now turn to a more quantitative analysis. An interesting feature 
of Bianchi-I spacetime lies in the comparison with the isotropic spacetime. It is important to note 
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FIG. 16: Variation of the number of e-foldings with the changing shear scalar at the bounce for > at 
the bounce. The values of cj) and a^ are in the units of rrip^ and Zpj respectively. 



that in the isotropic spacetime, if i;^(0) = 3.14r?T,pi at the bounce then the value of (f) at the bounce 
is fixed to |0(O)| ss 0.905 rripj as the bounce always occurs when energy density saturates its upper 
maximum p = Pmax = 0.41yOpi. Whereas, in the Bianchi-I spacetime, since the energy density may 
not be saturated at the bounce, |</>(0)| < 0.905m,p[ is allowed. In this way, there is a freedom in the 
choice of initial inflaton velocity at the bounce in Bianchi-I spacetime, as compared to the isotropic 
spacetime. Table-HIl shows the number of e-foldings for various initial values of the inflaton field with 
0(0) = (1.0, 2.0, 3.0, 3.14, 4.0, 5.0)?n,pi. For each of </>(0), as we just discussed, there are several 
possible initial <j) in the Bianchi-I spacetime, unlike in the isotropic spacetime where the inflaton 
velocity at the bounce is fixed for a given (j){0). This table shows the number of e-foldings for 
different initial velocities (/)(0) = (—0.002, 0.002, —0.02, 0.02)7TT-pj with varying initial shear scalar 
for each of the (/)(0) given at the bounce. The last two rows of Table -|TI]with |</>(0)| = 0.905 mpj 
correspond to the isotropic spacetime, characterized by o"^(0) = 0. It is evident from the analysis 
of Table-In] that, if the inflaton is initially rolling down, then, the Bianchi-I spacetime can produce 
significantly more number of e-foldings than in the isotropic spacetime. For example, if the initial 
value of inflaton field is (p{0) = 3.147Ti,pi then the number of e-foldings in Bianchi-I spacetime is 
of the order of 60 while for the same initial value of inflaton the isotropic spacetime produces 
e-foldings as small as 3.103, if the inflaton is initially rolling down. In the same way, the maximum 
number of e-foldings for 0(0) = S.Ompi in Bianchi-I spacetime is of the order of 50 compared to 
the isotropic spacetime where A^ ~ 1.648. In order to generate A^ ~ 60, in the isotropic spacetime 
for an inflaton initially rolling down, the initial value of the inflaton should be 0(0) ~ 5.50mpi 
whereas, the same amount of e-foldings can be generated in the Bianchi-I spacetime for a lower 
initial value of inflaton 0(0) < 5.50mpi (see Table-HH) . In this way, presence of anisotropy widens 
up the window of the values of the initial inflaton field, at the bounce, which produces a given 
amount of inflation. 



3. Phase portrait and attractor behavior 

Let us now explore the properties of the dynamical trajectories in LQC. We will discuss the 
trajectories of both the matter and the gravitational sector i.e. (0, 0) for the inflaton and (log(a) , H) 
for the gravitational part respectively. A phase portrait of Bianchi-I spacetime in LQC is shown in 
Fig.[T71 The left plot corresponds to the inflaton in which the dashed thick lines show the isotropic 
trajectories with the initial conditions 0(0) = ib3.147Ti-pi and 0(0) w ±0.905 TTipj at the bounce. The 
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1.0 
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3.0 


3.14 


4.0 


5.0 


-0.002 


0.002 
0.277 
0.967 
1.441 
3.184 
5.983 
6.831 


6.477 
6.486 
6.486 
6.487 
6.486 
6.485 
6.484 


23.548 
23.555 
23.559 
23.561 
23.562 
23.563 
23.564 


56.959 
56.971 
56.976 
56.979 
56.981 
56.983 
56.984 


62.404 
62.426 
62.428 
62.429 
62.428 
62.425 
62.424 


101.016 

101.03 

101.037 

101.041 

101.044 

101.046 

101.048 


155.575 
155.592 
155.601 
155.606 
155.61 
155.613 
155.615 


0.002 


0.002 
0.277 
0.967 
1.441 
3.184 
5.983 
6.831 


6.654 
6.645 
6.644 
6.644 
6.645 
6.646 
6.647 


23.887 
23.880 
23.877 
23.874 
23.873 
23.872 
23.871 


57.447 
57.436 
57.431 
57.427 
57.425 
57.424 
57.423 


62.833 
62.812 
62.809 
62.809 
62.811 
62.812 
62.813 


101.659 
101.645 
101.638 
101.633 
101.63 
101.628 
101.627 


156.365 
156.347 
156.339 
156.333 
156.33 
156.327 
156.325 


-0.02 


0.049 
0.287 
0.976 
1.450 
3.194 
5.993 
6.841 


5.765 
5.787 
5.795 
5.795 
5.789 
5.779 
5.776 


23.548 
23.555 
23.559 
23.561 
23.562 
23.563 
23.564 


56.959 
56.970 
56.975 
56.978 
56.981 
56.982 
56.983 


60.278 
60.352 
60.377 
60.381 
60.361 
60.349 
60.333 


101.016 

101.03 

101.037 

101.041 

101.044 

101.046 

101.048 


155.575 
155.592 
155.601 
155.606 
155.61 
155.613 
155.615 


0.02 


0.049 
0.287 
0.976 
1.450 
3.194 
5.993 
6.841 


7.409 
7.385 
7.376 
7.376 
7.382 
7.394 
7.397 


27.221 
27.233 
27.217 
27.203 
27.191 
27.182 
27.175 


59.486 
59.502 
59.479 
59.457 
59.439 
59.425 
59.415 


65.001 
64.924 
64.898 
64.893 
64.915 
64.924 
64.943 


105.226 
105.248 
105.218 
105.188 
105.164 
105.145 
105.131 


161.673 
161.700 
161.663 
161.626 
161.596 
161.572 
161.554 


-0.905 





12.595 


2.087 


1.648 


3.103 


17.997 


47.994 


0.905 





66.835 


111.65 


168.524 


177.454 


237.553 


318.805 



TABLE II: This table summarizes the variation of number of e-foldings with varying shear scalar at the 
bounce, for various initial (j) at the bounce. 

solid (blue) curve and the dotted-dashed (red) curves correspond to the LQC trajectories with the 
initial conditions at the bounce as 0(0) = ib3.14?TT-pi and 0(0) = ±0.002 rripj. In this phase portrait, 
the first and third quadrant correspond to an inflaton which is initially rolling up, while the second 
and fourth quadrant represent the trajectories when the inflaton is rolling down the potential. 
The slow-roll is characterized by the horizontal line, with an almost constant 0, to which all the 
trajectories meet in their future evolution. It is evident that for a rolling down inflaton with the 
given initial conditions, the Bianchi-I trajectory remains on the slow-roll curve for a longer period 
of evolution than the isotropic trajectory. In this way, the presence of anisotropic shear at the 
bounce (even a small value) leads the dynamical trajectory to meet the slow-roll curve sufficiently 
early which, in contrast to the isotropic case, results in higher number of e-foldings if the inflaton 
is initially rolling down the potential. 
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FIG. 17: The left figure shows the phase-space trajectory of the scalar field and the right plot is the plot of 
the trajectory in the gravitational phase-space (log(a) ,H). 



Analysis of the phase-space trajectories of gravitational variables log(a) vs. H gives information 
about the isotropization of the spacetime. The right plot in Fig.[T7]shows the comparison of isotropic 
trajectories with Bianchi-I trajectories for various values of initial shear at the bounce. The thick 
(black) curve shows the isotropic trajectory, the dashed (blue) curve corresponds to Bianchi-I 
trajectory with initial shear scalar o"^(0) = 5 x 10~^/pj , the dashed (orange) curve with initial 
shear (T^(O) = 5.7 x 10^^ /pj^ and the solid (red) line with (T^(O) = 6.81 /pj^ at the bounce, while 
the initial conditions on the inflaton are being fixed at 0(0) = — 0.002 ?7ipj and 0(0) = 3.14?7ipi. 
The black dot denotes the bounce point of the mean scale factor. It is clear from the plot that, 
close to the bounce, the evolution of Bianchi-I spacetime trajectories are different from that of the 
isotropic spacetime. Whereas, away from the bounce when the curvature is much smaller, then all 
the trajectories tend to meet the isotropic trajectory, leading the spacetime to isotropize. 

We now discuss, in detail, a feature of the attractor behavior of Bianchi-I spacetime in the 
effective description of LQC. Fig. [18] shows the evolution of trajectories starting with the same 
initial 0(0) = 3.14 mpi but different initial 0. For clarity, we only show the first and the fourth 
quadrant of the phase portrait. The initial values of for the trajectories in the fourth quadrant 
are 0(0) = (-0.905, -0.856, -0.727, -0.584, -0.536, -0.407, -0.289, -0.171, -0.053)m|i, and 
in the first quadrant, 0(0) = (0.0, 0.053, 0.171, 0.289, 0.407, 0.536, 0.584, 0.727, 0.856, 0.905)m2j. 
The left most trajectory, in the fourth quadrant, denoted by o"^ = corresponds to the isotropic 
spacetime with 0(0) = — 0.905 mpp and the right most trajectory in the first quadrant marked 
with (T^ = represents the evolution of isotropic spacetime with 0(0) = 0.905 m-pj. All the other 
trajectories correspond to Bianchi-I spacetime with non-zero initial anisotropic shear scalar. Recall 
that in the effective description of Bianchi-I spacetime the shear scalar has an upper bound, cr'^^^ = 
11.57 /p[ . Evolution of trajectories starting with initial a^ = (T^^ax ^^^ shown by solid (orange) 
curves denoted by "cr^ = 11.57". All the other Bianchi-I trajectories correspond to smaller values 
of initial shear scalar the bounce. The slow-roll curve is characterized by the horizontal line with 
an almost constant ~ 1.97 x lO^'^TTipj, to which all the trajectories with various initial 0(0) and 
(T'^(O) meet in their future evolution. Thus, like in the classical theory, the slow-roll curve turns out 
to be an attractor for all the trajectories shown in the phase portrait. This implies that irrespective 
of the initial anisotropic content of the spacetime, all the Bianchi-I spacetime trajectories in the 
effective description of LQC do meet the slow-roll curve in their future evolution, if the initial 0(0) 
is kept high enough. Thus, we conclude that the presence of non-vanishing anisotropic shear does 
not prevent the occurrence of slow-roll inflation with 0^ potential in Bianchi-I LQC. 
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FIG. 18: This figure shows the evolution of various trajectories starting with different 4>{Q) = 
(-0.905, -0.856, -0.727, -0.584, -0.536, -0.407, -0.289, -0.171, -0.053)to|,i starting in the fourth quad- 
rant and (/)(0) = (0.0, 0.053, 0.171, 0.289, 0.407, 0.536, 0.584, 0.727, 0.856, 0.905)m^i in the first quadrant 
(from left to right). The solid black curves marked with cr^ = correspond to the isotropic trajectories 
and the solid (orange) curves denoted by a^ — 11.57 (in Planck units) correspond to trajectories with max- 
imum shear scalar at the bounce. For all the trajectories, the initial value of the inflaton is kept fixed at 
0(0) = 3.14 mpi. All the values shown in this plot are in Planck units. 

In the effective description of LQC, the isotropic spacetime is characterized by two conditions: 
the first condition requires the vanishing of the shear scalar o"^ = at all times and the second 
condition requires the energy density to be p = pmax at the bounce. This leads to an interesting 
feature of effective dynamical trajectories of Bianchi-I spacetime regarding its isotropic limit. That 
is, in order to approach the isotropic limit of Bianchi-I spacetime in effective dynamics of LQC, 
the energy density at the bounce must tend to approach its maximum, i.e. p — )• Pmax (|';^(0)| — ^ 
0.905 rripj, if (/)(0) = 3.14mpi), along with the shear scalar vanishing, i.e. o"^(0) — )• at the bounce. 
This stands in contrast with the classical theory, where o"^ — )• is sufficient for Bianchi-I spacetime 
to approach the isotropic limit, irrespective of the value of the energy density. As discussed earlier, 
in the classical theory, for a fixed initial energy density, the isotropic limit can be approached by 
continuously varying the initial cr^ to zero as shown in Fig. [71 Thus, there is an isotropic trajectory 
corresponding to every initial energy density in the classical theory. In contrast, in LQC, the 
energy density for an isotropic spacetime must be p = Pmax at the bounce. In this way, for a 
given initial value of the inflaton field, there are several classical isotropic trajectories, but there 
are only two isotropic LQC trajectories (as shown in Fig. llSp . These two trajectories correspond 
to different signs of initial (j), such that p = pmax at the bounce. Therefore, for a given initial value 
of the infiaton field, there is a unique isotropic curve in each quadrant of the (j) — (p phase portrait 
of Bianchi-I spacetime in effective description of LQC. The two isotropic curves corresponding to 
(f)(0) = 3.14mpi are shown in Fig. 1181 one in first and the other in the fourth quadrant. Starting 
from any given trajectory in Fig. 1181 if one decreases (j){0) to — 0.905 mpj, which is the value of 4> 
at the bounce in the isotropic spacetime for an infiaton which is rolling down, then the Bianchi-I 
trajectories approach towards the isotropic trajectory in the fourth quadrant. While if one increases 
0(0) to 0.905 mpj, the value of 4> at the bounce in the isotropic spacetime for an infiaton which 
is rolling up, then the isotropic trajectory in the first quadrant is approached. In this way, the 
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FIG. 19: This plot shows the 3D phase portrait of various trajectories starting with different (/)(0) < in 
the near bounce region. The black dot marks the bounce point. The sohd (orange) line corresponds to 
the maximum shear scalar at the bounce a^ — cr^jax — 11.57 Zp?. The thick black line in the cr^ = plane 



corresponds to the isotropic trajectory. Clearly, as the value of initial 
trajectories approach the isotropic trajectory. 



tends to 



-0.905 mpj, the Bianchi-I 



isotropic spacetime turns out to be an "isotropic attractor" of Bianchi-I spacetime in the effective 
description of LQC. This behavior is easy to understand from a 3D plot of the phase trajectories 
as shown in Fig. ll9l where we show the evolution of Bianchi-I trajectory in the vicinity of the 
bounce. The solid black trajectory in the cr'^ = plane corresponds to the isotropic spacetime with 
0(0) = — 0.905 rripj and all the other ones correspond to Bianchi-I spacetime with different initial (j) 
and (T^ at the bounce. It is clear to see from this figure that as 0(0) — ?• — 0.905 mpj, the trajectories 
tend to approach the isotropic trajectory in the a^ = plane. 

It is important to note that the isotropic attractor behavior in LQC is obtained by varying 
both the initial energy density at the bounce and the shear scalar, whereas in the classical theory 
the isotropic spacetime can be approached for any fixed initial energy density. If one fixes the 
energy density, at the bounce in LQC, to any other value than /Omax by specifying 0(0) and 0(0) 
and tries to vary the shear scalar, then it turns out that the shear scalar can not be varied to a 
value less than a non-zero minimum value which depends on the energy density at the bounce. In 
this way, again, unless the energy density at the bounce tends to pmaxj Bianchi-I trajectories do 
not approach the isotropic spacetime. Thus, in LQC, isotropic attractor is not obtained for those 
Bianchi-I trajectories for which p ^ Pmax at the bounce. 

To obtain a complete set of initial data, for the trajectories in Fig. [181 we give Pi(0), ci(0), 0(0) 
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for varying (f)(0), and the values of C2(0), 03(0) are computed by solving the Hamiltonian constraint 
and the mean Hubble rate being, H = at the bounce. In the numerical simulations, it turns out 
that for a fixed value of ci(0), there is a fixed range of the values of (p{0) for which the Hamiltonian 
constraint is satisfied. That is, (j){0) can not be continuously varied from — 0.905 ?7ipj to 0.905 rripj 
for a fixed value of ci(0), in the numerical simulations performed. In this way, depending on 
various fixed values of ci(0), one can generate a continuous family of trajectories by varying (/>(0) 
within the corresponding range, which depends on the value ci(0). In each of these ranges, there 
will be a limiting curve with minimum shear scalar to which the other trajectories in the range 
tend, as the shear scalar in the initial data decreases. For the simulations shown in Fig. 1181 one 
has to choose atleast two different values of ci(0), one for |i;^(0)| < 0.5839 JTipj and another for 
\<j){0)\ > 0.5839 m|i. 

It is now clear from the discussion of the phase portraits of the infiationary Bianchi-I spacetime 
that: 

• Trajectories starting with a wide variety of initial conditions meet the slow-roll trajectory in 
the future evolution. This indicates slow-roll is an attractor for Bianchi-I spacetime. 

• Isotropic spacetime behaves like an "isotropic attractor" for Bianchi-I spacetime if the 
isotropic limit is taken by considering both cj^(O) — )• and p(0) — t- Pmax at the bounce. 
For the trajectories, for which the initial energy density does not approach pmax at the 
bounce, the isotropic limit does not exist. 

Hence, the two attractor behavior present in the classical Bianchi-I spacetime with cj)^ potential, 
is also present in the effective description of LQC, with the above noted subtlety in the way isotropic 
limit of Bianchi-I spacetime is taken. It is worth emphasizing that due to quantum geometric effects 
in the pre- inflationary era, the physical trajectories of the effective description are very different 
from those of the classical Bianchi-I spacetime. It is due to these differences that the conditions to 
approach the isotropic limit are modified in the effective description of LQC as compared to the 
classical theory. 

IV. DISCUSSION 

Though the inflationary paradigm has been successful in explaining the physics of the early 
universe, it is incomplete without understanding the way the initial singularity is resolved and the 
way inflation begins from generic initial conditions. In this work, we have analyzed these questions 
for the (j)"^ inflation in the setting of an anisotropic spacetime described by the Bianchi-I model 
using the effective spacetime approach of LQC and made a comparative analysis with the classical 
theory. The underlying quantum geometric effects originating from LQG result in a non-singular 
evolution. The backward evolution in LQC is devoid of the big bang singularity. The latter is 
replaced by a bounce of the mean scale factor when the spacetime curvature becomes Planckian. 
The energy density and the shear scalar remain bounded throughout the evolution in LQC. To 
our knowledge, this is the flrst such study of the inflationary scenario in Bianchi-I spacetime in a 
non-singular setting. 

One of the key questions in the study of the evolution of the universe is its isotropization from 
a generic anisotropic universe in the pre-inflationary era. Previous studies of Bianchi-I inflation in 
the classical theory and brane world cosmology established that the presence of anisotropy does 
not prevent inflation, and it is actually conducive to the number of e-foldings obtained during the 
inflationary phase. These studies for quadratic potential took the initial velocity of the inflaton to 
be negative (i.e. the inflaton initially rolling down the potential). It was found that the amount of 
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inflation for such initial conditions always increased with increasing shear in the classical theory 
as well as in brane- world scenarios [20|. We revisited the analysis in the classical theory by taking 
a more generic initial data by including the case when inflaton is initially rolling up the potential. 
In this case, we found that in comparison to the corresponding isotropic evolution, the amount 
of inflation actually decreases with an increasing anisotropic shear. This happens because of the 
greater anisotropy, the Hubble friction in the Klein-Gordon equation increases which turns up the 
decay of kinetic energy of the inflaton. Due to the faster decay of the kinetic energy, an inflaton 
which is initially rolling up, stops at a lower value of the potential which decreases the amount of 
e-foldings in the subsequent slow-roll evolution. It is worth emphasizing again that the presence 
of anisotropic shear does not rule out inflation. We showed that for accelerated expansion to 
take place, the potential energy must win over the kinetic and the shear energy. Owing to this 
property, in the inflationary era, shear term is always less than the potential energy. There may 
however be, a small period during the accelerated expansion when anisotropic shear contribution 
((T^/lGvrG) while being less than potential term {V{4>)), is still greater than the kinetic term {(j)'^/2). 
In this sense there may be a brief duration of shear dominance in the inflationary phase of Bianchi- 
I spacetime. We also showed by explicit simulations that if one starts with highly anisotropic 
initial conditions such that two of the scale factors are contracting and third one is expanding, 
in the inflationary phase all of the scale factors expand. This property can be attributed to the 
value of the anisotropy parameter ej which always decreases in an expanding universe. Therefore, 
irrespective of large initial ej, the shear scalar decreases and due to the accelerated expansion there 
is a point, in the future evolution, where ej < 2/\/3 which corresponds to all expanding directions 
of the Bianchi-I spacetime. Hence, isotropization is inevitable in an inflating classical Bianchi-I 
spacetime with a (j)"^ potential. 

In the analysis of the loop quantum effective dynamics of the inflating Bianchi-I spacetime, the 
quantum gravitational corrections come into play in the near bounce phase in the pre-inflationary 
era. Due to the presence of an upper bound on the anisotropic shear and the energy density, the 
set of initial conditions are quite different from those in the classical theory. Also, to give initial 
conditions on the bounce point, along with satisfying the Hamiltonian constraint, one has to make 
sure that the mean Hubble rate is zero, a condition which was not present in the classical theory. 
Our analysis shows that the amount of inflation in LQC does not change monotonically with the 
initial anisotropic shear at the bounce. It has a turning point at o"^ whose value depends on the 
initial value of inflaton and its velocity (but is independent of the sign of the initial velocity). For 
an inflaton initially rolling down the potential, the number of e-foldings increase with anisotropy 
only when initial o"^ is less than aj. On the other hand, if the inflaton is initially rolling up, then 
the number of e-foldings decrease with anisotropy only when initial o"^ is less than o"^. As in the 
classical theory, isotropization is inevitable in the inflationary LQC Bianchi-I spacetime starting 
from arbitrary anisotropic initial conditions. 

Analysis of the phase space trajectories demonstrate that there is a double attractor behavior in 
effective dynamics of LQC for Bianchi-I spacetime as in the classical theory, however one has to be 
careful in taking the isotropic limit at the bounce. In the classical theory, phase space trajectories 
approach the slow-roll attractor in the future evolution. This holds true also in LQC. There is also 
a non slow-roll isotropic attractor which is reached by taking the limit cr^ — >• for any value of 
energy density in the classical theory. In LQC, since energy density is saturated to its universal 
maximum value pmax at the bounce, the second attractor is reached only for those trajectories 
starting at the bounce which satisfy o"^ — )• along with p — )• pmax- 

We conclude our discussion by pointing out that in the Bianchi-I spacetime sufficient number 
of e-foldings can be obtained by starting at a significantly lower value of the inflaton field at the 
bounce than in the isotropic LQC when the inflaton is initially rolling down the potential. Further 
there is more freedom in specifying the initial data at the bounce in the anisotropic model. This can 
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be of significance in the study of observational signatures of quantum gravitational effects in the 
pre-inflationary stage of LQC. As an example, it has been recently suggested that for the isotropic 
LQC, a narrow window depending on the value of the inflaton field at the bounce exists for which 



there can be potentially distinct signatures of LQC in the cosmic microwave background 7% l73l|. 
Our analysis suggests that such a window potentially widens in the presence of anisotropies if the 
inflaton is initially rolling down in the pre-inflationary epoch. Owing to widening the window 
of the value of the inflaton field in Bianchi-I spacetime for the same number of e-foldings in the 
isotropic case, there could be potentially additional corrections which will again be of quantum 
gravitational in origin. A complete quantum treatment of the field perturbations in Bianchi-I 
anisotropic background is expected to bring out additional possible observational signatures of 
loop quantum effects in the primordial power spectrum. 
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